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We investigate the ballistic motion of electrons in III-V semiconductor quantum wells with Rashba 
spin-orbit coupling in a perpendicular magnetic field. Taking into account the full quantum dynamics 
of the problem, we explore the modifications of classical cyclotron orbits due to spin-orbit interaction. 
As a result, for electron energies comparable with the cyclotron energy the dynamics are particularly 
rich and not adequately described by semiclassical approximations. Our study is complementary to 
previous semiclassical approaches concentrating on the regime of weaker fields. 
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I. INTRODUCTION 

The coupling between the orbital and the spin degrees 
of freedom of itinerant carriers in semiconductors is a ma- 
jor direction of work in today's spintronics research. An 
early key example is given by the proposal of a spin field- 
effect transistor put forward by Datta and Das already 
in 1990i. The merit of this paradigmatic theoretical con- 
cept is that it allows, once realized, all-electrical control 
of electron spins in two-dimensional III-V semiconductor 
structures, avoiding any magnetic field. On the other 
hand, at about the same van Houten et al. published 
a pioneering both theoretical and experimental study of 
magnetic focusing of electrons in semiconductor quan- 
tum wells^. Here the control over the orbital degree of 
freedom of the carriers is achieved by a perpendicular 
magnetic field of typically moderate strength, very anal- 
ogously to the classical cyclotron motion. This effect has 
been demonstrated both for conduction-band electrons 2 - 
and for valence-band holes 3 -^; very recently cyclotron or- 
bits in an electron focusing experiment were directly im- 
aged using scanning probe microscopy^. 

Other recent experimental studies have addressed the 
question whether it is feasible to spatially separate (and, 
in turn, separately detect) carriers in different spin-spilt 
subbands of a quantum well via such transverse focusing 
tcchniques&i^. In these investigations, the spin split- 
ting of subbands was either provided by a strong in-planc 
component of the magnetic ficld^ , or, more relevant for 
the present study, by specific contributions to spin-orbit 
coupling acting on the spin of carriers in semiconduc- 
tor quantum wells^^. In the latter case, different initial 
spin states provide via spin-orbit interaction a separation 
of carriers in real space, an effect which has already at- 
tracted also significant theoretical interes t 10 ' 11 : 12 ' 13 : 14 ' 15 . 
In the present paper we provide a complementary theo- 
retical study of cyclotron motion and magnetic focusing 
under the influence of spin-orbit interaction using a fully 
quantum mechanical approach. 

A related phenomenon is the predicted zitterhewegung 
of carrier wave packets in the presence of spin-orbit 
couplingi^^d^^^^^^^. Zitterhewegung 



of free electrons described by the four-component Dirac 
equation was originally predicted by Schrodinger and oc- 
curs for wave packets which contain solutions of the free 
Dirac equation of both positive and negative energy^. 
In effective models for electrons and holes in semicon- 
ductors the intrinsic spin-dependent energy splitting due 
to spin-orbit coupling can lead to a similar oscillatory 
zitterhewegun g 16 ! 11 'i 21 ' 25 . However, the latter effect is pre- 
dicted to occur on time and length scales being much 
more favorable for experimental detection compared to 
the situation of free electrons where zitterhewegung has 
never been observed so fa r 16 ' 17 . On the other hand, a uni- 
fying aspect of these two phenomena is given by the fact 
that zitterhewegung of itinerant band carriers in semicon- 
ductors occurs due to spin-orbit interaction which can be 
viewed as the nonrelativistic of the strong coupling be- 
tween spin and momentum being manifest in the Dirac 
equation. 

Moreover, the interplay between spin-orbit coupling 
and cyclotron motion in a perpendicular magnetic field 
was already studied theoretically in some detail in 
Refs . 25 ' 27 . Here the authors concentrate on semiclassical 
approximations, and on an analogy between the Jaynes- 
Cummings model of atomic transitions in a radiation field 
and the Rashba Hamiltoniar>22. in a perpendicular mag- 
netic field, an aspect to be briefly reviewed below. In 
the present paper we report on numerical evaluations of 
the full quantum mechanical dynamics of a free electron 
in a two-dimensional quantum well with spin-orbit inter- 
action and a perpendicular magnetic field, avoiding any 
further approximation. As explained in the appendix, 
our approach is so far technically limited to the Hilbcrt 
space of the first few ten lowest Landau levels. For typi- 
cal electron energies of a few meV, this restriction corre- 
sponds for usual III-V semiconductor materials to mag- 
netic fields of a few tesla. Such fields are somewhat larger 
than those considered in circumstances of semiclassical 
approximations neglecting Landau quantization, and in 
this sense our present study is complementary to those 
previous investigations. For definiteness we will also con- 
centrate on spin-orbit coupling of the Rashba type, al- 
though also other effective coupling terms can be consid- 
ered. Finally we note that a complementary theoretical 
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study of conduction-band electrons being subject to spin- 
orbit coupling and a homogeneous in-plane electric field 
was given very recently in Ref.— . 

This paper is organized as follows. In section |TT] we 
summarize the essential properties of the Rashba model 
in a perpendicular magnetic field. We discuss the anal- 
ogy to the Jaynes Cummings model of quantum optics, 
and we describe in detail the initial states used for the 
numerical simulations of time evolutions to be discussed 
in section HITl All further technical details can be found 
in the appendix. We close with conclusions in section [TVl 



II. MODEL AND APPROACH 



connect different Landau levels. Note that the Hamil- 
tonian (including the spin-orbit part) can be expressed 
in terms of a and a + only, no further orbital operators 
occur. Therefore its eigenstates have the same Landau 
level degeneracy as in the absence of spin-orbit coupling. 

Fixing a certain intra-Landau-lcvel quantum number, 
we denote by |n, a) = {{a + ) n / \/n\) |0, a) a state in the 
n-th Landau level with spin direction a G {T,4}- Then 
|0, t) is an eigenstate with energy e = (hw c + <7/zb-B)/2, 
and all other eigenstates are of the forr n 30 i 31 
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with energy 



We consider an electron in an n-doped quantum well 
being subject to Rashba spin-orbit coupling 3 ^ and a ho- 
mogeneous perpendicular magnetic field coupling both to 
the orbital degrees of freedom as well as to the spin, i.e. 
the single-particle Hamiltonian reads 
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Here m the the effective band mass, tt = p + eA/c is 
the two-component kinetic momentum with the canoni- 
cal momentum p and the vector potential A generating 
the magnetic field B along the growth direction of the 
quantum well chosen as the z-axis, B = Vxi. The ef- 
fective Rashba spin-orbit coupling parameter is denoted 
by a, g is the effective g-factor, [i B the Bohr magneton, 
and a are the usual Pauli matrices. Note that the Rashba 
Hamiltonian can be viewed as a momentum-dependent 
field coupling to the electron spin, an interpretation we 
will use later on in the discussion of numerical results. 
Moreover, in the following we will assume, without loss 
of generality, that the product of the electron charge 
(— e) = — |e| and the magnetic field strength B is always 
positive, (—e)B > 0, i.e. B points along the negative 
^-direction. 



A. Spectrum and eigenstates 

Defining the usual bosonic operators 
1 I 



fulfilling [a, a + ] = 1 and £ = ^Jhc/\eB\ being the mag- 
netic length, the Hamiltonian reads 
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where uj c = \eB\/(mc) is the cyclotron frequency, and we 
have defined — o x ± ia y . The operators a and a + 
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£„ = hu c n ± y 2n—^-hw c + - (hjj c + gfisB) (5) 
and the amplitudes parametrizing the eigenstates read 
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Thus, the energy levels and eigenstates of the system 
are characterized by the interplay of three energy scales: 
The cyclotron energy e c = huj c , the Zeeman energy sz = 
gfisB, and the Rashba energy er = ma 2 /h 2 . 



B. Analogy to the Jaynes-Cummings model 

As it was recognized recently in Ref^, the Hamilto- 
nian ^ is formally equivalent to the Jaynes-Cummings 
model for atomic transitions in a radiation field. This 
model has been studied very intensively in theoretical 
quantum optics, and the time evolution of orbital and 
spin operators has been obtained in terms of analytical 
but rather implicit expression s] 32 ' 33 . To explore this anal- 
ogy it is useful to separate the Hamiltonian into two com- 
muting parts, H = Hi + H.2, with 



Hi = fv^i c ( a + a + 



1 +CT 2 



(8) 
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where mo is the bare electron mass. Then the time evolu- 
tion of the position operators in the Hciscnbcrg picture, 



f H (t) = e iHt/h f ^ e -im/h 



(10) 
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can be written a a 25 ' 32 i 33 

XH(t) + WH{t) = x + iyo 
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where the operator-valued frequencies ui± are given by 



hoj± 



-Hz± 



ma" . 
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(12) 



and xo, yo ar e the usual coordinates of the center of the 
classical cyclotron orbit which commute with the Hamil- 
tonian and are therefore constant in time. 

The result (fTTj) is correct but still not very explicit. In 
particular, the operator character of the quantities uj± 
poses severe obstacles against evaluating this expression 
for a given initial state. Therefore, in the present work 
we follow a different route towards the full quantum dy- 
namics by expanding the initial state of the system in 
terms of its eigenstates. 



C. Gauge and initial state 

In the very general considerations so far there was not 
any necessity to specify the gauge of the vector poten- 
tial A. For the practical calculations to be described 
below, however, we shall work in the Landau gauge 
A = (0, Bx, 0) where the spinless orbital eigenstates in 
the absence of spin-orbit coupling have the following well- 
known form 
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labelled by a wave number k corresponding to transla- 
tional invariance in the y-direction, or, equivalently, by a 
guiding center coordinate kl 2 for the x-direction. H n (x) 
are the usual Hermite polynomials, and the phases of 
the above wave functions have been adjusted to fulfill 
a\n, k) = y/ri\n — 1, k). 

In what follows we will be interested in the quantum 
dynamics of an initial state \rp) being a direct product of 
an orbital and a spin state, 



(14) 



where the spinor components are related to the usual 
polar angles d, ip of the initial spin direction via K = 
exp(-«V/2) cos(i?/2), A = exp(«V/2)sin(i?/2). As a 
generic initial orbital state we consider 
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i.e. a normalized Gaussian wave packet of spatial width d 
and initial momentum ftfco along the y-axis, i.e. the direc- 
tion of translational invariance of the Hamiltonian. The 
initial position of the particle is at the origin, (ip\r\ip) = 0. 
The energy of the above initial state can be expressed 

as 
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where the bar denotes complex conjugation. Conversely, 
for a wave packet of the above form with fixed energy 
E and given width and spin state, the initial momentum 
reads 
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D. Time evolution 

A conceptually straightforward way to evaluate time- 
dependent expectation values is to expand the initial 
state in terms of the eigenstates of the above system and 
use the matrix elements of the desired operator in this 
eigenbasis. For instance, for the kinetic-momentum op- 
erator, this approach formally reads 
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where we have already anticipated that these operators 
are diagonal in the intra-Landau-lcvcl quantum number 
k; the same holds for the spin operators a nit). The 
summation over /i;, i G {1, 2} runs over fii = ± for n, > 
and fa =| for ?i; = 0. From the expectation values of 
the position operators can be obtained from those of the 
momenta via 
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(19) 
(20) 



where the expectation values of the constant centers of 
the classical cyclotron motion are given by (ip\xo\ijj} = 
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ko£ 2 , (ip\yo\ij}} = 0. As explained in detail in the ap- 
pendix, the integration over the wave numbers k can be 
performed separately and serves as an input for the nu- 
merical evaluation of the remaining sums. For any fur- 
ther technical details, we refer the reader to the appendix. 



III. RESULTS 

We now present the results of numerical simulations 
of the time evolution of expectation values described in 
Eq. (UHJ). All relevant technical details can be found in 
the appendix. In all simulation we assume the Rashba 
coefficient to be positive, a > 0. 



A. Cyclotron motion 

Let us first investigate the influence of spin-orbit cou- 
pling on the cyclotron motion in general. Fig. [1] shows 
the the particle orbit evaluated in terms of the expecta- 
tion values (tp\rH(t)\tp) =■ (rff(i)) of a wave packet of 
initial width d = 1. 01 and group wave number ko = 2.0 jl 
for various initial spin states. The Rashba energy is 
eji = 0.2hiu c while the Zccman energy is, for simplicity, 
put to zero here. In the left (right) top panel, the spin 
points initially along the positive (negative) x-direction. 
The middle and bottom panels show the corresponding 
data for the y- and z-direction, respectively. The total 
simulation time is always t = 30/w c . The strictly circular 
motion (dotted lines) with radius kgi 2 occurring in the 
absence of spin-orbit coupling is shown in all graphs as a 
guide to the eye. The magnetic length £ can conveniently 
converted into practical units via t = 257 Aj ^B/Tesla. 

All six graphs have the appearance of a more or less 
distorted spiral. The prima vista most regular mo- 
tion is found in the two top panels where the initial 
spin direction is collincar with in initial direction of the 
momentum-dependent coupling to the spin described by 
the Rashba Hamiltonian. This situation was investi- 
gated very recently in Rcf. 27 in the framework of several 
schemes of semiclassical approximations. In one of these 
approaches the spin is assumed to follow in an adiabatic 
fashion the momentum-dependent field coupling to it, 
where both quantities are taken to be classical variables. 
We will discuss below to what extend this approximation 
leads to useful results in the parameter regime considered 
here where the energy of the initial wave packet is com- 
parable with the cyclotron energy. 

In Fig. [2] we have plotted the corresponding spin dy- 
namics expressed in terms of the time-dependent expec- 
tation values ((?#(£)) . The initial conditions in the three 
panels are the same as in the left column of Fig. [1] The 
solid lines show the modulus of the vector |(erjy(£))| of 
the time-dependent expectation values of spin compo- 
nents. This quantity can be used as a measure of entan- 
glement between the electron spin and its orbital degrees 
of freedo m 34 ! 35 . In fact, when tracing out the real-space 



degrees of freedom, the time-dependent reduced density 
matrix of the spin reads 
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(21) 



with eigenvalues X± (t) = (1± | (<?#(£)) |)/2. Thus, a mod- 
ulus of |(<7ij(i))| = 1 (as present in the intitial condition 
at t = 0) corresponds to a direct product of spin and or- 
bital state with a reduced spin density matrix of rank 1, 
while a vanishing modulus | ((?#(£)) | = indicates max- 
imal entanglement between spin and orbital degrees of 
freedom, and the reduced spin density matrix is propor- 
tional to the unit matri x 34 ! 35 . As seen in Fig.[2J the mod- 
ulus | (<7H{t)) | is generically clearly smaller than unity sig- 
nalling between spin and real-space coordinates, an effect 
certainly beyond semiclassical approximations. 

Let us now come back to the investigation of the adi- 
abatic semiclassical approximation employed in Rcf. 2 '. 
Here both spin and particle momentum are treated as 
classical variables, and the projection of the spin on 
the instantaneous direction of the momentum-dependent 
field is assumed to be constant. Intuitively, this assump- 
tion corresponds to strong spin-orbit coupling 2 ^, as it is 
the case for a Rashba energy of er = 0.2huj c studied 
above. To investigate the validity of this adiabatic ap- 
proximation, we introduce 



Qi(t) 

and 



\(*H(t))\ 



(22) 



02W: - \(*«{t))WH(t))\ ^ 

The first quantity is the projection of the vector (<?#-(£)) 
onto the direction of the momentum-dependent field eval- 
uated in terms of (??#(£)), whereas in Q 2 we have addi- 
tional divided by \(&H(t))\ in order to eliminate the ef- 
fects of entanglement discussed above. For the adiabatic- 
scmiclassical approximation to be valid, Q\ and Q 2 
should be reasonably constant in time. 

Fig. [3] shows the time dependence of Qi and Q 2 for 
the same system parameters as in Fig. [TJ In particu- 
lar (tth(0)} points along the positive y-direction which 
means that the momentum-dependent field coupling to 
the spin is initially in the x-direction in spin space. In 
the top panel, the spin points initially along the posi- 
tive (negative) x-direction with Qi(0) = Q 2 (0) = +1 
(Qi(0) = Q2(0) = —1). The middle and bottom panels 
show the analogous data with the spin initially aligned 
along the y- and z-axis, respectively. Here we have al- 
ways Qi{0) = Q2(0) = 0. As seen in the figure, Qi{t) 
and Q2{t) significantly deviate from a constant value even 
in the case where the spin is initially fully aligned with 
the momentum-dependent field (top panel). From these 
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observations we conclude that this adiabatic scmiclassi- 
cal approximation is rather problematic in the parameter 
regime studied here where the total energy of the electron 
wave packet is of the order of the cyclotron energy. In 
fact, the behavior of the system is only rather poorly rep- 
resented by introducing two different cyclotron radii cor- 
responding to two spin directions as suggested in Refill. 
On the contrary, the dynamics are much richer and show 
trajectories reminiscent of chaotic behavior. The lat- 
ter observation becomes even more significant at smaller 
initial group wave number fc as one can see in Fig. 0] 
where we have plotted trajectories analogous to those in 
Fig. [T] but with a shorter initial group wave vector of 
only fco = 0.5 /£. For lower cyclotron energies, however, 
we expect the semiclassical approximation of Refill to 
be significantly better fulfilled than in the regime studied 
here. The smaller the cyclotron energy compared to the 
total energy of the wave packet, the larger the number of 
Landau levels to be included in the numerical simulation. 
As explained in the appendix, such simulations require 
the precise numerical evaluation of high-order Hermite 
polynomials, a task which technically limits our approach 
to the regime where the cyclotron energy is comparable 
with the energy of the initial wave packet. In this sense, 
our study is complementary to previous semiclassical ap- 
proaches concentrating on the regime of weaker fields. 

Finally Fig.[5]shows the orbital dynamics for again the 
same system as in Fig. [T] for various values of the initial 
group wave number fco an d the spin initially always point- 
ing along the positive x-axis. For a better comparison the 
the components of (r) are given in units of k £ 2 . Clearly, 
the dynamics become more "regular" (or less "chaotic" ) 
the larger the initial group wave number fco. 

In the above investigations we have concentrated on 
the position operator rjj (t) to describe the time evolution 
of the initial state chosen as a Gaussian wave packet. 
Regarding the width of this wave packet (as opposed to 
its center (Fij(t))) let us consider the case of vanishing 
spin-orbit coupling. Here the components of the position 
operator are straightforwardly obtained as 

x H (t) = x(0) + -^-sh-i(uj c t) 
rnujc 

+ (i-cosM)) , (24) 

muj c 

y H {t) = 2,(0) + ^ sin (u c t) 
mu) c 

_7Tx(0) , l _ coa , Uc ty ( 25 ) 
muj c 

(26) 

and are completely analogous to the classical cyclotron 
motion. For an initial state Gaussian wave packet given 
in Eq. (|15[) the time dependent width reads 

(x%(t))-(x H (t)) 2 + {y%(t))-(y H (t)) 2 
= d 2 + \ - d 2 ) (1 - cos (w c t)) . (27) 



Thus, differently from the dispersive dynamics of a wave 
packet in the absence of a magnetic field, the width does 
not increase to infinity but remains bounded and rather 
oscillates with the cyclotron frequency, similarly to the 
time evolution of a coherent state in a harmonic oscilla- 
tor. In the presence of spin-orbit coupling we expect the 
time evolution of the width of the initial state to be more 
complex but still essentially bounded. 



B. Magnetic focusing 

Let us now turn to the issue of magnetic focusing under 
the influence of spin-orbit coupling. A magnetic focus- 
ing experiment is conceptually very simple and sketched 
in FigO Electrons enter a quantum well at a location 
x = Xi, follow a ballistic cyclotron orbit, and impinge 
again on the boundary of the system at a location x = Xf. 
In the absence of spin-orbit coupling, this difference in co- 
ordinate depends, just as in the classical case, only on the 
applied magnetic field and the initial group wave number 
fco, Xf — Xi = 2k £ 2 . Thus, using an appropriately located 
detection contact one can study electron transport as a 
function of these two quantities. 

In a typical experiment, however, rather the energy E 
of electrons (defined by the Fermi energy of the injecting 
lead) not their momentum is fixed with both quantities 
being connected via Eq. (|I7[) . Therefore, the wave num- 
ber fco of an injected electron will in general depend on 
its spin state; only in the absence of both Zeeman cou- 
pling and spin-orbit interaction fco is independent of the 
electron spin. In turn, for random initial spin directions, 
Xf — Xi will be distributed according to some probability 
density P x {xf — Xi). FigJT] shows a numerical evaluation 
of P x (xj — Xi) for a wave packet of width d = £, total and 
energy E = 2.0Hlo at zero Zeeman coupling and different 
Rashba energies. The initial angular coordinates cosd 
and ip determining the complex amplitudes k and A in 
the initial state (|14p were chosen at random from uniform 
distributions in the interval [—1,1] and [0, 2ir], respec- 
tively. The data is averaged over 500000 randomly cho- 
sen initial spin states each. At a small Rashba energy of 
only er = m(a/h) 2 = O.Olftw (top left panel) the dynam- 
ics depend only very weakly on spin and P x (xf — Xi) is 
strongly peaked around Xf — Xi = 2k £ 2 « 3.16, the clas- 
sical cyclotron diameter expected in the absence of spin- 
orbit coupling. With increasing Rashba energy, this peak 
undergoes a broadening with the maximum of the prob- 
ability density being located at smaller values of xf — x%. 
Fig. [5] shows the same type of data at a Zeeman energy 
of ez — O.lhuj which does not lead to any qualitative 
difference. 

In summary, the initial narrow peak of the probability 
density P x (xf — xi) at small spin-orbit coupling broad- 
ens with increasing Rashba energy and develops an non- 
trivial structure in terms of a maximum at small argu- 
ments with a broad shoulder reaching to higher values. 
The structures seen in Figs. 17181 appear to be somewhat 
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different to the results of Ref^i where a splitting of the 
conductance peak as a function of magnetic field was 
found for increasing Rashba coupling. These two peaks 
can be related to two different effective cyclotron radii 
corresponding to two initial spin states with respect to a 
quantization axis being perpendicular to the initial mo- 
mentum and the magnetic field, an observation similar 
to the "strong coupling" semiclassical scenario of Ref^i. 
Thus, in the light of these investigations one could also 
expect a double-peak structure to develop in the proba- 
bility density P x (xf — Xj). However, the investigations of 
Refi^i work at clearly higher electron energies compared 
to the cyclotron energy, and Landau quantization is ex- 
plicitly neglected. This is in contrast to the present study 
which works at larger cyclotron energies taking into ac- 
count the full quantum dynamics of the problem. 



IV. CONCLUSIONS AND OUTLOOK 



employed^. It is an interesting question whether the in- 
clusion of both kinds of couplings leads to qualitativly 
new observations. A particular situation is reached if 
both terms occur with the same magnitude, where, for 
zero Zeeman coupling, a new conserved spin operator 



arise; 



,36 



Further possible extensions of the present work in- 
clude the study of valence-band holes (as opposed to 
conduction-band electrons) with an effective spin-orbit 
coupling being trilinear in the momentum, and electron 
or hole dynamics under the influence of an additional 
in-plane electric field. 
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We have studied the ballistic motion of electrons in 
III-V semiconductor quantum wells with Rashba spin- 
orbit coupling and a perpendicular magnetic field. Dif- 
ferently from previous investigations, our numerical ap- 
proach takes into account the full quantum mechanics 
of the problem and is technically limited to situations 
where the cyclotron energy is of the same order as the 
energy of the initial electron wave packet. For typical 
experimental parameters, this restriction corresponds to 
magnetic fields of a few tesla. Such fields are larger than 
those considered previously in circumstances of semiclas- 
sics neglecting Landau quantization, and in this sense our 
present study is complementary to such semiclassical ap- 
proaches. As a result, in the parameter regime considered 
here the electron dynamics are particularly rich and not 
adequately described by semiclassical approximations. 
An interesting issue for futher investigations here in- 
cludes the question whether the seemingly "chaotic" tra- 
jectories shown in section [til Al are truely ergodic. More- 
over, it is tempting to attribute the irregularity of these 
trajectories to the zitterbewegung predicted previously for 
electron motion in two-dimensional electron gases with- 
out magnetic & e i^SJLiSJS^L2^2^2S£US.. What 
both phenomena have indeed in common is the fact that 
they are the result of spin-orbit coupling, and the irregu- 
lar motion of electrons in a perpendicular magnetic field 
is the consequence of the non-equidistant spectrum of 
Landau levels induced by spin-orbit interaction. 

In this study, we have concentrated on spin-orbit cou- 
pling of the Rashba type. However, the situation of linear 
Dresselhaus coupling can be treated analogously^ as it 
only couples, like the Rashba term, pairs of neighbor- 
ing Landau levels with opposite spin. If both types of 
spin-orbit coupling terms are present, all Landau levels 
are coupled, and the single-particle Hamiltonian cannot 
be diagonalizcd analytically anymore. In this case the 
eigensystem of the Hamiltonian needs to be computed 
numerically, or appropriate approximations have to be 



APPENDIX A: TECHNICAL DETAILS 

1. Overlap with basis states 

The overlap of the initial orbital state (fTS)) and the ba- 
sis states (flU]) of the usual Landau levels can be expressed 
as 



(n,fc|0) = H)" 



x cxp 




2 \ n/2 



(Al) 



(A2) 



where [x] denotes the largest integer not larger than x, 
and the second of the above equations shows explicitly 
that the overlap is well-behaved at d = £. The above 
expressions can be obtained by using the explicit form of 
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the Hcrmitc polynomials, 



d" 



H n {x) = (-l)»e» 



E 

p=0 



{-If 



p\ (n — 2p)\ 



(2x) 



n — 2p 



(A3) 
(A4) 



Finally, the overlap of the initial state (IM|) with the spin- 
ful eigenstates Q is given by 

<0, /c, T 1-0) = K(0,k\<f>), (A5) 
{n,k,n\ijj) = u%K(n,k\<t>) +v%\(n-l,k\<f>) . (A6) 



The in-plane components of the spin operator have the 
matrix elements 



(0, k, T \a x \n,k,fx) = uf5i, n , 
(0, fc, t |cr a |n, fc,/i) = -iuf(5i. n , 



(All) 
(A12) 



2. Matrix elements 

As already stated, the matrix elements of the kinetic 
momentum as well as the spin operators are diagonal 
with respect to the wave number k. For the kinetic mo- 
mentum, the matrix elements read explicitly 



(0,k,1 \% x \n,k,fi) = j-j=u^5i^ n , (A7) 



(0,fe,t \ix y \n,k,ii) 



(ni, k, Hi\TT x \n 2 , k,fJ, 2 ) 

= (V^7TT« +1 + V^«+i) 



h 

(ni,k,fJ,i\Tr y \n2,k,fi2) 

— ih 

X IV5 ni '" 2_1 



ih 



>— S n 1 ,n 2 + l 



(A8) 



(A9) 



(A10) 



(ni,fc,/zi|cr x |n 2 ,£;,^2) 



l u n 1 +l u ni,n 2 



+ «£<+i*ni,n a +l , (A13) 
(ni,fc, i m|a- !/ |n2,fc,M2) = -ra£^ +1 <J ni ,n 2 -i 

whereas a z is diagonal in the Landau level index n, and 
the nonvanishing matrix elements read 



(o,fc,r k 2 |o,fc,t) = i, 

(n, k,^i\a z \n, fc,/x 2 ) 



(A15) 
(A16) 



3. Explicit time evolution 



Using the expressions given in the previous sections, 
the time-evolved expectation values of the components 
of the kinetic momentum can be formulated as 



8 



((7t x )h(*)> = Rc J y/2- £ etto-'V* (|«| 2 |<| a Jo + *A«I ) 



M^E E 



n=l /ii,/i2=± 



|«| 2 (-v/n+TK'HWn+ll 2 + ^^n^r+lCl) J « 



+|A| 2 (v^+I^< i < 2 +1 <^ + V^km'K+iI 2 ) ■/«-! 
+sa (Vn+TK 1 | 2 <;id + idi 2 ) 4 



(A17) 



((%)*(*)) = Re i ^ e^"^ 4 (| K | 2 |<| 2 Jo + kA«/ ) 



Re. 



-^?E E e * 



n— 1 /ii ,/j,2— ± 



| K | 2 (^Am|< 1 | 2 |< 2 + ll 2 + MCCiCi) J « 



+SA (^Kf^!^! + V^Whdl 2 ) In 



(A18) 



where we have defined 

/oc 
dk((f)\n,k)(n,k\(j)) , 
-oo 

/•OO 

J„ = / dk{(j)\n, k){n + 1, , 

dk{(f>\n- l,k)(n+ l,k\<f>) 



Using Eq. (|A2|) it is straightforward to derive explicit 
expressions for these integrals in terms of finite sums to 
(A19) be evaluated numerically. As an example, for /„ one finds 

(A20) 
(A21) 



In = 



I 2 



d 2 

7? p,q=0 



f—iy+q n i2"+ 1 - 2 (p+9) 
p\q\ (n- 2p)\{n-2q)\ 



d 



2 \ P+1 



\ 2n-(p+q) 
-jT ] Rn-(p+q)(£,d,ko) 



(A22) 



with 



R m (e,d,k ) = -= \ 



\2 



dk{kl) 2m exp I - (J^pp + d2<yk ~ fc °) 2 



eX P - ,4 , ,2,2 ,4 -7 , ,2,2 , M dk ° -( A2 3) 



^ + d 2 P + d 4 J v \ i 4 + d 2 e 2 + d 4 J V 4/ zm \\ e 4 + d 2 i 2 + d 4 

1 

Note that, despite the imaginary argument of the Her- Thus computing the quantities /„ ( and similarly J n and 
mite polynomiali^m, Rm is always real and positive. 
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K n ) requires again the evaluation of Hermitc polynomi- 
als. At large Landau level index n this operation limits 
the accuracy of the present numerical approach. Given 
the data /„, J n , K n , the summations (|A17[) . (|A18[) are to 
be performed numerically where the sum over the Lan- 
dau level index n can be truncated at a sufficiently large 
energy. For the simulations presented in this work it is 
sufficient to take into account the first 25 Landau levels 
where the evaluation of Hermite polynomials is numeri- 
cally unproblematic. 

It is noteworthy that the quantities /„, J„, K n fulfill 
certain sum rules which provide a convenient check on 
numerical evaluations. For instance, normalization of the 
initial state \4>) obviously requires 

oc 

7 « = 1 • ( A24 ) 

n=0 

I 



+ 



<K(*)> = 



+ 



<(*h(*)> = 
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n=0 
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Analogously one derives 

00 1 r ]2 i 

£ yfifc+7jK n =— a¥ ~ - 2 k\e . (A26) 

71 = 

Finally, coming back to physical expectation values, 
the time-evolved spin components read 
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FIG. 1: Orbital dynamics of a wave packet of initial width 
d — l.Oi and group wave number fco = 2.0/1 for various ini- 
tial spin states. The Rashba energy is er = 0.2ftoj c while the 
Zeeman energy is put to zero. In the left (right) top panel, the 
spin points initially along the positive (negative) s-direction. 
The middle and bottom panels show the corresponding data 
for the y- and z-direction, respectively. The simulation time 
is always t = 30 /u c . The stricly circular motion (dotted lines) 
with radius kol 2 occurring in the absence of spin-orbit cou- 
pling is always shown as a guide to the eye. 
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FIG. 2: Spin dynamics as expressed in terms of the time- 
dependent expectation values corresponding to the left col- 
umn of Fig. [1] The solid lines show the quantity |(<th(£))| 
which is a measure of entanglement between the spin and the 
orbital degrees of freedom. 
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FIG. 3: The quantities Qi (solid lines) and Q2 (dashed lines) 
defined in the text as a function of time for the same sys- 
tem parameters as in Fig. [1] In the top panel, the spin 
points initially along the positive (negative) s-direction with 
Qi(0) = Qa(0) = +1 (Qi(0) = Q 2 (0) = -1). The middle and 
bottom panel show the analogous data with the spin initially 
aligned along the y- and z-axis, respectively. Here we have 
always Qi(0) = Q 2 (0) = 0. 
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FIG. 4: Orbital dynamics for the same system as in Fig.[T]but 
with a smaller initial group wave number of only ko — 0.5/1. 
Again, the stricly circular motion occurring in the absence of 
spin-orbit coupling is always shown as a guide to the eye. 
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FIG. 5: Orbital dynamics for the same system as in Fig. Q]for 
various values of the initial group wave number ko and the 
spin initially always pointing along the positive x-axis. For 
a better comparison the the components of (r) are given in 
units of kol 2 . 
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FIG. 6: Schematic sketch of a magnetic focusing experiment. 
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FIG. 7: The probability density P x {Xf —Xi) for a wave packet 
of width d = £, total and energy E = 2.0/ioj at zero Zeeman 
coupling and different Rashba energies. The data is averaged 
over 500000 randomly chosen initial spin states each. Note the 
different scale of the y-axis in the top, middle, and bottom 
panels. 
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FIG. 8: The probability density Px(xf — Xi) for the same 
situation as in Fig[7]but with a Zeeman energy of e~ = O.lficj. 



